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RESUME 

Ce  rapport  d€crit  un  modele  mathSmatique  pour  la  solution  du 
probleme  de  propagation  de  faisceaux  laser  se  d6pla£ant  dans  1' atmosphere 
turbulente  et  corriges  par  une  optique  adaptable.  Celle-ci  est  simulee 
&  l'aide  d'une  relation  mathematique  simple  mais  suffisamment  generale 
pour  repr6senter  la  plupart  des  systemes  existants.  Le  modele  tient 
compte  egalement  des  effets  de  vibration  du  faisceau  et  de  coherence 
partielle  de  la  source.  La  methode  permet  de  pr§dire  les  profils  de 
l'intensite  moyenne  et  de  la  variance  de  1' intensity  pour  des  niveaux 
arbitraires  de  scintillation.  On  presente  quelques  solutions  typiques 
pour  des  faisceaux  de  longueur  d'onde  de  3.8  et  10.6  ym.  Ces  resultats 
illustrent  les  performances,  apres  propagation  dans  la  turbulence,  d'un 
systeme  d'optique  adaptable  en  fonction  de  la  longueur  d'onde,  du 
nombre  d' elements  actifs  et  du  niveau  de  bruit.  (NC) 


ABSTRACT 

This  report  describes  a  mathematical  model  for  solving  the 
propagation  problem  of  laser  beams  traveling  in  atmospheric  turbulence 
and  corrected  by  adaptive  optics.  The  modeling  of  the  adaptive  optics 
is  mathematically  simple  but  sufficiently  general  to  encompass  the 
majority  of  the  existing  systems.  The  model  also  includes  the  effects 
of  beam  jitter  and  partial  coherence  of  the  source.  The  method 
allows  the  prediction  of  the  average  irradiance  and  the  irradiance 
variance  beam  profiles  for  arbitrary  scintillation  levels.  Typical 
solutions  are  presented  for  3.8-and  10.6-ym  laser  beams.  These  results 
illustrate  the  performance,  after  propagation  in  turbulence,  of  an 
adaptive  system  as  a  function  of  the  wavelength,  the  number  of  active 
elements  and  the  noise  level.  (U) 
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1.0  INTRODUCTION 

Atmospheric  turbulence  seriously  affects  the  propagation  of 
laser  beams .  The  refractive  index  turbulence  induces  random  phase  and 
amplitude  modulations  which  cause  the  beam  to  scintillate,  to  lose  its 
spatial  and  temporal  coherence,  to  wander  about  its  axis,  and  to  spread 
out.  A  prediction  model  for  the  average  irradiance  and  the  irradiance 
variance  profiles  of  laser  beams  traveling  in  such  a  turbulent  medium 
was  developed  at  the  Defence  Research  Establishment  Valcartier  (DREV) 
and  is  described  in  Refs.  1,  2  and  3.  The  principal  advantages  of  the 
model  are  that  it  is  applicable  at  arbitrary  scintillation  levels  and 
that  the  solutions  can  be  computed  by  straightforward  finite-difference 
techniques.  The  results  of  Refs.  1  to  3  show  that  the  numerical  solutions 
agree  very  well  with  the  experimental  data  and  that  the  classical  per¬ 
turbation  solutions  are  recovered  in  the  weak  scintillation  limit. 

In  this  report,  we  modify  the  DREV  model  to  include  the  corrective 
effect  of  adaptive  optics  which  has  now  been  proven  quite  efficient 
in  compensating  for  wave  front  distortions  induced  by  turbulence. 

The  many  aspects  of  this  new  technology  are  well  reviewed  in  Ref.  4. 

The  technique  will  likely  have  an  important  impact  on  several  optical 
and  infrared  applications  in  the  atmosphere.  It  is  therefore  essential 
for  the  design  and  analysis  of  upcoming  systems  to  be  able  to  simulate 
adaptive  optics.  The  existing  numerical  codes  in  the  literature 
concern  mainly  the  simulation  of  the  optical  systems;  propagation  is 
mostly  reduced  to  computer-generated  random  input  signals.  This 
report  addresses  itself  to  the  intrinsic  propagation  problem.  The 
constitutive  closure  relations  of  Refs.  1  to  3  are  rederived  to  model 
the  operation  of  adaptive  optics .  The  approach  neglects  technical 
particularities  but  this  is  done  for  the  sake  of  generality  and  no 
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problem  is  foreseen  for  dealing  with  specific  systems.  The  purpose  is 
to  demonstrate  the  efficiency  of  the  DREV  model  to  simulate  the 
propagation  of  adaptively  corrected  laser  beams  in  turbulence. 

Section  2.0  reviews  the  basic  features  of  our  propagation  model. 
Section  3.0  describes  the  modifications  made  to  simulate  adaptive 
optics,  beam  jitter,  and  partial  coherence  of  the  source.  Section  4.0 
proposes  a  finite-difference  algorithm  and  Section  5.0  presents  and 
compares  calculation  results  at  two  wavelengths,  A  =  10.6  pm  and  3.8  pm. 

This  work  was  performed  at  DREV  between  October  1979  and  June 
1980  under  PCN  33B07,  Atmospheric  Propagation  of  Laser  Beams. 


2.0  PROPAGATION  MODEL 


Our  mathematical  model  for  beam  propagation  in  turbulence  is 
described  in  Refs.  1,  2  and  3.  The  model  is  in  the  form  of  a  closed 
set  of  simultaneous  partial  differential  equations  for  the  first-  and 
second-  order  statistical  moments  of  the  complex  wave  amplitude.  From 
the  solution  of  these  equations,  the  average  irradiance  and  the  irradiance 
variance  beam  profiles  can  be  easily  calculated. 

The  basic  features  of  the  model  are  recalled  from  Ref.  3.  First, 
the  instantaneous  random  scalar  electric  field  E  of  the  assumed  mono¬ 
chromatic  wave  propagating  in  the  z -direction  under  negligible  polariza¬ 
tion  effects  is  written  as  follows: 


E  =  A  exp  [jk(z  +  <(>)  -  jwt]. 


[1] 
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where  A  is  a  complex  amplitude,  k$  is  a  phase  function,  k  =  nQu/c  is 
the  optical  wave  number,  nQ  is  the  unperturbed  index  of  refraction, 
to  is  the  optical  angular  frequency  of  the  source,  c  is  the  speed  of 
light  in  free  space  and  j  =  V^T.  Second,  eq.  1  is  substituted  for  E 
in  the  scalar-wave  equation  which,  under  the  paraxial  approximation 
and  after  separation  for  A  and  $,  becomes 


( 

( 


V 

V 


vj  V  =  V(N  -  n0)/n0  , 

V)  A  +  jA  V  ■  V-^|-v2A  =  0, 


[2] 

[3] 


where  N  is  the  instantaneous  random  index  of  refraction,  J7  is  the 
gradient  operator  in  the  plane  normal  to  the  o-z  axis,  and  V  =  ,24>.  In 
obtaining  eq.  2,  it  was  assumed  that  |N-n  |  /nQ«l,  which  is  well 
verified  in  atmospheric  turbulence.  Third,  the  random  functions  are 
written  as  sums  of  an  average  and  a  fluctuating  part,  i.e. 


N  =  <N>+n;<n>=0, 
V  =  <V>+v;<v>=0, 
A  =  <A>  +  a;<a>=0, 


[4] 

[5] 

[6] 


where  the  pointed  brackets  denote  ensemble  averaging.  Then,  the  equa¬ 
tions  for  the  first-  and  second-  order  statistical  amplitude  moments 
are  easily  derived  from  eqs.  2  to  6,  but  they  contain  more  unknowns 
than  there  are  equations.  This  is  the  classical  closure  problem  which 
affects  the  treatment  of  turbulent  phenomena  governed  by  nonlinear  or 
quasi-linear  stochastic  equations  such  as  eqs.  2  and  3.  For  the  present 
application,  this  problem  was  solved  in  Ref.  3  where  the  unknown  moments 
are  mathematically  related  to  the  first-  and  second-order  amplitude 
moments,  thus  closing  the  set  of  equations.  For  a  wave  originally 
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focused  at  z  =  F  in  a  transversely  statistically  homogeneous  medium, 
i.e.  £  <  N  >  =  V,  <  v.v  >  =  0,  the  resulting  moment  equations  are 

±<A->-  +  —L^  •  v  <  A  >  +  >  K(Z)<A>-RU):W<A>-4rv^<A>=0, 

3z  z  —  F  z  —  F  2  2k 

[7] 

—  aa*  —  +  — r—  •  v  <  aa*  >  +  <  -  *  >- -lReal/R(z)l  :  w  <  aa*  > 

3z  z-F  z  —  F  2  IJ 

yF—  v2<aa*>=  2Real{(R(z)- v<  A>)- v<  A>*| 

+  Real  | K(z) l  <  A  >  <  A  >  *  ,  [8] 

-  <  aa  >-  +  — i—  •  v  <  aa  >  +  ------  +  j  T(z)  <  aa  >  — !-  R(z) :  VV  <  aa  > 

dz  z  —  F  z  —  F  2 

— 1_  vi  <  aa  >  =  2(R(z)  -  V<A>)-v<A>  +  K(z)  <  A  >  2 .  [9] 


where  0  is  the  diffractional  far-field  half-angle  beam  divergence.  For 
a  Gaussian  beam  with  an  irradiance  e -folding  radius  rQ,  0  =  l/krQ.  The 
coefficients  J&(z)  ,  K(z)  and  T(z),  which  constitute  the  basis  of  the 
closure  model ,  are  given  by 


R(Z)  =  ^rJ ^Tj^8  <  v(z-r)v(u.s)  >  exp 


—  2u  |  F  —  u 
F  —  z 


r  —  s 


O  —  oo 

Z  co 


2 

|] 


[10] 


:<v(Z,r)v<u,s)>exp[-2^  Fy_u2u|  r~s  |] 


T(z) 


-{ 


-5  Z<ZA- 

F  -  z  z 


[11] 


r,  F  -  z/2  z 

C  — - -  ;  z  >  zA 

F  -  z  za2 


[12] 
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where  C  =  0.15  -  0.01  j  is  a  numerical  empirical  constant  and 
zA  =  n  is  the  propagation  scale  characteristic  of 

turbulence  fading.  It  is  noted  that  the  semi -empirical  eq.  12  is 
not  required  to  solve  for  the  average  irradiance  but  only  for  the 
irradiance  variance.  This  empirical  input  could  be  omitted  only  at  the 
expense  of  a  much  higher  level  of  complexity,  namely  the  solution  for  the 
amplitude  covariance  function  which  satisfies  a  partial  differential 
equation  of  at  least  two  more  independent  variables.  Hence,  in  view  of 
the  simplicity  and  accuracy  of  eq.  12,  as  demonstrated  by  the  results 
of  Refs.  1  to  3,  the  need  for  the  universal  empirical  constant  C  is 
considered  acceptable. 

The  quantities  of  interest  are  the  average  irradiance 


<I>  =  <A><A>*  +  <  aa*  >, 


[13] 


and  the  irradiance  variance  which,  under  the  approximation  of  Gaussian 
statistics  for  the  complex  amplitude  A  (e.g.  Ref.  5),  is  given  by 

or  =  <  aa*  >-  +  <  aa  >  <  aa  >*  +  2  <  A  >  <  A  >*  <  aa*  > 

[14] 

+  <A><A><aa>*  +  <A>*<A>»<aa>. 


The  coefficients  g,(z)  and  K(z)  are  functions  of  the  covariance 
of  the  random  vector  v.  The  stochastic  equation  for  %  is  easily  derived 
from  eqs.  2  to  5.  For  an  originally  spherical  phase  front  focused  at 
z  =  F,  it  is  given  by 


+  — - —  •  vv  +  — —=?  +  v'7v  =  Vn/n0  . 
dz  z  -  F  z  -  F 


[15] 
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which  is  the  paraxial  form  of  the  eikonal  equation  of  geometrical  optics. 
Hence  v(z,£)  is  the  component,  in  the  transverse  plane,  of  the  unit 
vector  parallel  to  the  geometrical  ray  passing  through  the  point  (z,  r) 
or,  in  other  words,  v(z,  r)  is  the  vector  angle  subtended  by  the 
geometrical  ray  at  (z,  jr)  . 

In  Ref.  3,  the  expressions  for  jl(z)  and  K(z)  were  derived  from 
the  solution  of  eq.  15,  with  the  boundary  condition  ^(o,  jr)  =  0.  In 
the  present  report,  we  will  rederive  these  expressions  to  account  for 
the  action  of  an  adaptive  optical  system  that  produces  random  fluctuations 
of  v  at  z  =  0  in  response  to  some  corrective  criterion.  Therefore,  our 
task  will  be,  first,  to  determine  the  statistics  of  ,v(o,  r) ,  second, 
to  solve  eq.  15  with  this  new  boundary  condition,  and  finally,  to 
derive  the  resulting  expressions  for  R(z)  and  K(z).  The  solution  to 
the  propagation  problem  is  then  straightforward  with  eqs.  7  to  9  forming 
a  closed  set  of  equations  easily  solved  by  finite  difference  techniques. 


3.0  MODELING  OF  R(z)  AND  K(z) 


In  this  section,  we  derive  the  coefficient  R(z)  and  K(z)  for  an 
adaptively  corrected  beam  propagating  in  homogeneous  Kolmogorov's 
turbulence.  Also  derived  are  the  contributions  to  JR(z)  and  K(z) 
resulting  from  partial  coherence  of  the  source  and  residual  jitter  of 
the  transmitting  optics. 

3.1  Adaptive  Optics  Model 

The  adaptive  optical  system  is  assumed  to  derive  its  information 
from  the  radiation  returned  by  an  unresolved  glint  on  the  target 
illuminated  by  the  outgoing  beam  at  focal  distance  z  =  F.  The  geometry 
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of  the  returned  signal  is  therefore  that  of  an  originally  spherical  wave 
and  is  schematically  illustrated  in  Fig.  1.  It  is  evident  that  an  auxil¬ 
iary  point  source  could  also  replace  the  reflecting  glint. 

The  quantity  of  interest,  or  the  information  available  to  the 
adaptive  system,  is  the  state  of  distortion,  at  the  transmitter  plane, 
of  the  phase  of  the  returned  spherical  wave.  In  the  notation  used  here, 
this  is  measured  by  the  random  ray  angle  v  which  satisfies  eq.  15. 

Since  |v|«l  for  atmospheric  turbulence,  the  nonlinear  term  in  eq.  15 
can  be  neglected  to  give 


FIGURE  1  -  Schematic  representation  of  the  reflected  wave  from  a 
glint  on  the  target 
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where  the  variables  £  r'  are  defined  in  Fig.  1.  Solving  with  the 
boundary  condition  v(o,  r'J  =  0,  we  find  in  the  transmitter  plane 

F 

»■("-{  dfKvn(f'r'= I  o  • 

o 

The  principle  of  operation  of  an  ideal  adaptive  correction  system 
is  as  follows.  A  sensor  measures  in  some  form  or  other  the  random 
phase  front  angle  v^(r),  a  processor  then  reverses  the  sign  of  vt(r) 
and  feeds  the  information  to  a  phasor  array  that  continuously  changes 
the  phase  front  of  the  outgoing  laser  beam  to  make  ^(o,r)  =  - v, (r) .  If 
the  latter  condition  is  realized,  the  reversibility  of  the  ray  paths 
predicts  that  the  phase  front  of  the  outgoing  beam  will  gradually 
become  undistorted  and  spherical  as  z  ■*  F.  Hence,  the  irradiance  profile 
of  the  corrected  beam  should  approach  the  diffraction  limit. 


In  practice,  the  phase  sensor  and  corrector  have  a  discrete 
number  of  active  elements.  Therefore,  the  system  cannot  reproduce  the 
function  -^(r)  on  every  point  of  the  aperture  of  the  outgoing  beams 
and  some  form  of  spatial  averaging  must  be  taken  into  account.  The  true 
relation  between  the  returned  function  v^fz)  and  the  theoretical  eq.  17 
is  particular  to  each  system  but  it  is  believed  that  the  following 
linear  expression  is  sufficiently  general  to  accommodate  a  wide  class 
of  applications: 


d2p  W,(p)  {vb(r,  +  p)  -  H(P)  Vt(r,  +  p) } . 


[18] 


where  r^  is  the  position  vector  of  the  i^  element,  v^(x^)  is  a  random 
function  representing  the  noise  of  the  system,  I\(g)  is  the  correlation 
coefficient  between  the  true  and  the  returned  wave  front  angle  for  the 
i^  active  element,  W^(p)  is  the  window  or  weighting  function  associated 
with  the  it^1  element,  and  a.  is  its  surface. 

’  l 
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The  solution  of  the  linearized  eq.  IS  for  the  outgoing  phase 
front  with  boundary  condition  1/(0 ,  r)  =  v,  (jr)  is 


o 


Using  eqs.  17  to  19  and  noting  that  £  =  F  -  g,  we  have 

vJJ d,‘,w'«"v''(Tb  '  fiS) 

<JX 

d2"  w^>  vn  (f-T=7 r>  +  trLp) 

<T,  O 

z 

* 4fdf7rivn(i'-f^7'')  '  l“] 

O 


An  expression  for  the  covariance  of  v  is  easily  derived  from  eq. 

20.  The  method  is  well  documented  in  Ref.  3.  The  simplifying  hypo¬ 
theses  are  that  the  correlation  length  of  the  turbulent  refractive 
index  is  much  smaller  than  the  propagation  distances  of  interest  and 
that  the  noise  v^  is  uncorrelated  with  the  index  n.  Both  these  approxi¬ 
mations  are  well  justified  in  practice.  The  resulting  solution  for 
Kolmogorov's  turbulence  is  then 


<  v(z.  r, )  v(u 


■5>,>  =  ,f-zhf--u)  iM-  W,<p' W'm  <  v"(  7777 r  ”  p)  '••>(  f~Ts  *  r)  > 

<7,  (Tl 

[f-ViTp -  tt.  '  -  'V  '<■  ’] 


dfiF-n1  w  rr-f  F-f 


F-  u 


^JdVw,(p)i»J’-Fz)(F_u) 

17,  O 

. . »-/.■ 


C  dflF-  n1  „  p  -  s'  F  -  S'  1 
J  (F  —  zMF  —  u)  M  LF-z  r  F  —  u  s  J 


[21] 


. . .  • 
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C  is  the  index- structure  constant,  i0  is  the  inner  scale  of  turbulence, 
is  the  unit  dyad  (or,  in  tensor  notation,  the  two-dimensional  Kronecker 

'N/ 

deltc)  and  a  is  the  unit  vector  along  £. 


The  covariance  of  as  defined  by  eq.  21,  is  a  discrete  function 
of  the  position  vectors  r^  and  of  the  active  elements  labeled  i  and 
j.  This  is  very  inconvenient  to  handle.  Thus,  to  extend  further  the 
analytical  treatment,  we  render  eq.  21  continuous  by  assuming  that  the 
domain  of  the  spatial  integration  over  cn  (or  a.)  remains  symmetrical 
about  any  point  r  (or  £) .  In  other  words,  we  assume  that  the  averaging 
effect  of  the  active  elements  is  invariant  with  position  on  the  aperture. 
Finally,  to  simplify  the  results,  the  integration  is  performed  over  a 
circular  region  of  effective  radius  rg  given  by 

re  =(S/mir)y\  [23] 


where  S  is  the  total  area  of  the  aperture  and  m,  the  number  of  active 
elements.  This  approximation  appears  reasonable  for  systems  with 
continuous  phase  correctors  such  as  deformable  mirrors,  but  it  should 
also  give  acceptable  results  even  for  segmented  mirrors  especially  when 
the  number  of  active  elements  is  large. 

As  shown  in  Ref.  3,  eq.  22  for  £(£)  is  an  empirical  expression 
connecting  two  theoretical  asymptotic  formulas  valid  respectively  for 
p  «  i.  and  p  »  SLq  .  It  is  therefore  only  consistent  to  proceed  in  the 
same  fashion  to  determine  the  spatial  averages  of  £(£) :  first,  compute 
the  averages  in  the  two  asymptotic  limits  and  then,  connect  the  results 
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by  a  single  empirical  continuous  expression.  For  the  present  application, 
the  weight  and  correlation  functions  W^(p)  and  I\(p)  are  taken  to  be 
unity.  Defining 


F-f  F-f 

A  - - —  r - —  s  , 

F - z  F-u 


[24] 


F  -  r 

we  have  for  A  >>  - —  r 


M[  +  l: 


M**'  = 


if)j5  d2r  m 

a 

[a  +  F  (p-T) 

1  2.43  C2 

J 

1 

A 

f  0-p  M 

a 

2.43  C2 
"  n2 

1 

A1/. 

[«-jaAaA] 

ld!™ 

] 

2.43  C2 
< 

1 

A'/. 

[6“7a^] 

]  •  [25] 


[27] 


and  for  A  << 


F  -  5 


e’ 


M 


d2pJT dV{  + 


2.43 Cn  (f-d2 

1  n  1  a1 


[28] 


Performing  the  integrations  in  eq.  28  and  using  o  =  irr^, 


M,'-1—  2,4n3,C"  5  |(.35f0)2  +■  0. 199 (F  r 
Similarly,  we  have 


[29] 


}  |  (.35f0)2  +  0.335 (FF;f)2r2|~,> 


[30] 
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Then,  the  two  asymptotic  regions  are  connected  by  a  single  algebraic 
formula  as  follows: 


M(A,L,)  = 


2-43  C2 
n« 


\ 

[a’  +  L,2]* 


3 


A2 

A2  +  L2, 


[31] 


where  i  =  1,  2,  3  and 

L2  =  (.350*  +  0-»99  (F-~/)2  t l ,  [32a] 

L2  =  L12=(.35f«)J  +  0.335  ^-=j^-r2  •  [32b, c] 

Finally,  from  eqs.  21  and  31,  the  covariance  of  v  is  given  by 


<  v(z,  r)  v(u,  s)  >  =  (-F— f(p  l-u)  dJpJ  d2r  <  vb (  +  p )  v^^s  +  r)  > 

F 


G  G 


*  J  M  <&  1 L'»  -  J  “w-tv-i, M  L’> 

o  o 


u 

J 


u  u 

*J  dfrp-ViF -u)HlA-‘->  •  [«) 


where,  for  symmetry  of  notation,  we  have  defined 

L2  =  (.35f0)2 


[32d] 


The  contribution  due  to  noise  depends  on  the  covariance  of  which 
is  characteristic  of  the  particular  system  under  study.  For  the  present 
general  application,  it  is  convenient  to  assume  a  homogeneous  and  iso¬ 
tropic  noise,  Gaussianly  correlated  with  a  e-folding  radius  equal  to  rg, 
i.e. 


^Jd!pJd!T<V,'(F^rr  +  P)  V”(f^T,"T)  >=  *2  4'XP{  “I 


F-  i 
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The  parameter  is  the  noise  strength,  i.e.  the  spatially  averaged 
variance  of  the  random  angular  error  of  the  active  elements.  Equation 
34  simplifies  the  algebra  and  should  give  results  representative  of  a 
real  system. 

The  expressions  for  the  coefficients  R(z)  and  K(z)  can  now  be 
obtained  by  substituting  eqs.  33  and  34  in  eqs.  10  and  11.  The  integration 
over  is  performed  as  described  in  Ref.  3.  For  the  noise  contributions, 
the  complete  integrals  are  done  analytically.  The  results  are  expressed 
as  follows: 


R(z)  =  Rb(z)  +  R,(z)  -  Rj(z)  -  R3(z)  +  R4(z) .  [35] 

K(z)  =  K  b  (z)  +  K.  ,(z)  -  K  2(z)  -  K  ,(z)  +  K4(z) ,  [  36  ] 


where  ^(z)  and  K^(z)  are  algebraic  expressions  of  elementary  tran¬ 
scendental  functions  to  be  given  later  in  nondimensional  variables  and 
where 


Ri(z) 


Ki(z) 


2-43  C-5 
"o 


3-24  Cjj 

nl 


t  Yi 


4d„(V ,  (F~  Q1  »'•>,) 

J  J  W’(F  -z)(F  +  z- 

o  o 


2u) 


(%) _ 

z  -  2u)  ' 


[37] 


[38] 


with 


Yi  =  F  ;  Y2  = 

=  z  ;  Y,  =  V4  =  u  . 

[39] 

—  jk  L2j  {F  —  u)(F  +  z  —  2u) 

[40] 

T'  2(z  —  u)  (F-f)2 

G(y)  =  y- 

1  +  3y  +  [  5/12  -  y  -  3r  ]  ^  r<5/6.  y)  . 

[41] 

X 

2 

II 

1  +  [  5/6  -  7  ]  ^F(5/6,  y)  , 

[42] 
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r(5/6,y)  is  the  incomplete  Gamma  function  defined  by  eq.  6.5.20  of 
Ref.  6. 


The  coefficients  R(z)  and  K(z)  of  eqs.  35  and  36  constitute  the 
principal  result  of  this  report.  They  permit  the  closure  of  the 
amplitude -moment  equations  and  they  model  the  effects  of  turbulence 
and  of  adaptive  optics. 

3.2  Jitter  Model 


Another  phenomenon  that  contributes  to  beam  divergence  is  jitter. 
Jitter  arises  from  random  vibrations  of  the  transmitter  platform  which 
cause  a  random  tilt  of  the  outgoing  phase  front.  We  consider  here  the 
residual  jitter  which  is  left  uncorrected  by  the  adaptive  system. 

To  determine  the  contributions  of  jitter  to  £(z)  and  K(z),  we 
proceed  in  a  manner  different  to  that  used  in  the  preceding  section. 

The  reason  is  that  the  phase  and  amplitude  fluctuations  resulting  from 
jitter  are  correlated  over  the  complete  field  of  propagation  which 
invalidates  a  basic  hypothesis  leading  to  eqs.  10  and  11.  We  therefore 
return  to  the  defining  constitutive  relations  given  in  Ref.  3,  i.e. 


<  aV'v  >  = 

-  K(z)  <  A  >, 

[43] 

<  av  >  = 

-  R(z):  V  <  A  >, 

[44] 

<  aav  >  = 

-  ^R(z):  V  <  aa  >, 

[45] 

<  aa*v  >  = 

-  ^  Real  {  R(z)  }:  V  <  aa*  >, 

[46] 

from  which  we  will  derive  j|(z)  and  K(z)  by  substituting  the  specialized 
jitter  solutions  for  %  and  a. 
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Let  a  be  the  tilt  angle  due  to  jitter  measured  from  the  axis  of 
propagation.  In  the  notation  of  this  report,  we  thus  have  in  the 
transmitter  plane 

v„(0.r)  =  tte„,  [47] 


where  the  subscript  a  refers  to  the  jitter  contribution  and  is  a 
unit  vector  of  random  direction  in  the  plane  normal  to  the  propagation 
axis.  Then,  from  eq.  19,  the  solution  ^  (z ,  x)  is 


v„(z.  r)  = 


F-z 


««„  • 


[48] 


Equation  48  shows  that^fz,  £)  is  uniform  in  the  transverse  plane, 
i.e.  independent  of  the  position  r. 

From  Ref.  3,  the  implicit  solution  for  the  fluctuating  amplitude 
caused  by  jitter  is 

*" u' rl ” irf T^irjf d’’ y" s) ' Vi Alu' !)”‘l { • 

°  ■“  (49] 

Since  •  v^  =  0,  the  amplitude  fluctuations  are  caused  by  the  random 

angular  motion  of  the  beam  as  a  whole.  In  other  words,  jitter  produces 
amplitude  modulations  on  the  same  scale  as  that  of  the  average  beam. 
Hence,  we  can  assume  that  the  scale  of  A,  at  least  the  part  correlated 

with  v  ,  is  greater  than  vT^-u)  X,  so  that  V  k{u,s)  can  be  taken  out  of 

^  F  -u  . 

the  integral  over  s  and  evaluated  along  the  line  £  =  -  r,  i.e. 

F-z 


V, 


A(u,  s)  =  — — “  Vt  A 
F  -  u 


(, 


F  —  u 
F-z 


[50] 


This,  in  fact,  is  the  geometrical  optics  approximation.  In  the  frame¬ 
work  of  this  approximation. 
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A 


F-u 

F-z 


0 


F-  z 
F-u 


A(z,  r) , 


[SI] 


and  eq.  50  becomes 


Vi  A(u.  s)  = 


(F  -  z)2 
(F-u)2 


VrA(z.r)  • 


[52] 


Finally,  using  eqs.  48  and  52  in  eq.  49  and  performing  the  integration, 
we  find 


a„  (z,  r)  =  -  a  z  e„  •  V  A(z,  r)  • 


[53] 


Equations  48  and  53  are  the  jitter  solutions  for  g  and  a  which 
can  now  be  used  to  derive  an  expression  for  <^a>a-  Multiplying  eq.  48 
by  53  and  taking  the  ensemble  average,  we  find 


<  va  >  „  = 


-  Fz 
(F-z) 


<  «2  e,^,,  •  VA  > 


[54] 


Assuming  that  is  statistically  independent  of  a  and  A  as  is  the  case 
in  most  applications  and  neglecting  third-order  moments  between  a  and 
A  which  are  0(<a3>),  we  can  simplify  eq.  54  as  follows: 


<  va  >„  = 


-  Fz 
(F-z) 


<  «2  >  <  e„e(t  >  •  v  <  A  > 


[55] 


Finally,  if  is  uniformly  distributed,  ce^  ga>  =  5/2  and  eq.  55 
becomes 


-  Fz  <  a2  > 


"  (F-z)  2 


6  •  v<  A  >  ■ 


<  va  > 


[56] 


W  1 
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To  derive  the  expression  for  <vaa>a>  we  multiply  eq.  53  by 
a^(z,r)  and  take  the  ensemble  average.  Using  eq.  48  for  anc* 
assuming  that  a  and  A  are  near  Gaussian  random  functions,  we  find  to 
the  same  order  of  precision  as  for  eq.  56 

<  vaa  >„  =  <  “  —  g  •  V<  aa  >  •  t57] 


Similarly, 


<  vaa*  >„  = 


-  Fz 
(F-z) 


<  Q2  > 

4 


6  •  v<  aa*  >, 


[58] 


and,  since  3,'X^  *  °» 


<  av  ■  v  >(t  =  0 . 


[59] 


Therefore,  from  comparison  between  eqs.  43  to  46  and  50  to  59,  it 
follows  that 


R  (z)  5 . 

<,U'  (F-z)  2 

[60] 

K„(z)  =  0  • 

[61] 

The  coefficients  R  (z)  and  K  (z)  model  the  effect  of  beam  jitter. 
«a  a 

The  parameter  =  <a2>  is  the  variance  of  the  uncorrected  tilt  angle 
of  the  outgoing  phase  front  caused  by  transmitter  vibrations. 

3.3  Partial  Coherence  Model 


Real  laser  sources,  especially  high-power  lasers,  do  not  generally 
have  complete  coherence.  They  are  partially  coherent  which  causes  a 
beam  divergence  greater  than  the  diffraction-limited  divergence.  As 


UNCLASSIFIED 

18 


shown  in  Ref.  7,  the  property  of  partial  coherence  is  modeled  by  a 
statistically  homogeneous  random  phase  over  the  source  plane.  With 
the  notation  of  eq.  1,  this  is  mathematically  represented  as  follows: 

A(0,r)  =  i/KrTexp  {j  'l'c(r) }  ,  [62] 


where  I(r)  is  the  irradiance  profile  of  the  source.  The  degree  of 
coherence  is  measured  by  the  correlation  function  of  A(o,r),  i.e. 

<  A(0,  r)  A«(0,_«)>_=  <  {j  (q,c(r)  _  vfrc^))}  >  •  [63] 

[I(r)  I(s>]* 

If  we  choose  that  V  has  a  Gaussian  probability  density  function  with 
zero  mean,  as  is  done  in  Ref.  7,  eq.  63  becomes 

<  A(O.r)  A*(M_>_  =  exp  /  -  J-<  (*c  (r)  -  *c(s))’  >  \  •  [  64 ] 

[l(r)l(s)p  l  2  J 

It  is  also  customary,  e.g.  Refs.  7-10,  to  consider  a  Gaussian  correlation 
function  for  A(o,r),  i.e. 

<  A(0, r)  A*(0,  s)  >  [  |r_s|2'[ 

[I(r)  I(s )]*  =CXP\  ~W~S  '  [65] 

where  r£  is  called  the  coherence  radius  of  the  source.  Although  the 
hypotheses  of  Gaussian  statistics  for  ¥  and  of  a  Gaussian  correlation 
function  for  A  are  not  essential,  they  permit  a  simplified  analytic 
treatment  of  the  phenomenon.  Hence,  from  eqs.  64  and  65  it  follows  that 


<  (*c  (r)  -  'I'c(s))2  >  =  • 


[66] 
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The  random  phase  H1  of  eq.  62  can  be  considered  a  boundary 
condition  to  the  geometrical  phase  kd>  of  eq.  1.  Since  the  equations 
of  the  present  model  deal  with  the  gradient  of  the  phase,  we  define,  in 
accordance  with  our  previous  notation,  the  random  vector 


vc(0,  r)  = 


rr  V  ♦c  (r)  • 

k 


[67] 


Then,  from  eq.  66,  we  obtain 


<  vc(0,r)vc(0,s)  > 


[68] 


We  note  from  eq.  68  that  the  covariance  of  x.c(o,t)  is  uniform  in  the 

transverse  plane,  i.e.  independent  of  £.  This  situation  is  mathemati- 

2  2 

cally  identical  to  jitter.  Therefore,  by  replacing  <a2>  by  2/k  r£  as 
indicated  by  eqs.  47  and  68,  the  results  of  the  preceding  section  are 
applicable  without  further  modifications.  Thus,  with  reference  to 
eqs. 60  and  61,  the  effect  of  partial  coherence  is  modeled  by  the 
coefficients 


RCU)  = 


Fz 

(F-z) 


[69] 


Kc(z)  =  0. 


[70] 


The  degree  of  coherence  of  a  source  is  often  defined  as  the 
ratio  q  of  its  far-field  divergence  to  that  of  a  fully  coherent  source 
with  the  same  irradiance  profile.  It  can  easily  be  shown,  for  a  Gaussian 
source  with  an  irradiance  e-folding  radius  equal  to  rQ,  that 


[71] 
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This  simple  relation  provides  a  practical  means  of  estimating  the 

coherence  radius  r  of  a  real  source, 
c 


4.0  METHOD  OF  SOLUTION 

We  have  derived  in  the  preceding  sections  the  relations  required 
to  close  the  system  of  partial  differential  equations  for  the  solution  of 
the  average  irradiance  and  the  irradiance  variance  of  optical  or  infra¬ 
red  beams  propagating  in  turbulence.  The  model  includes  the  effect  of 
adaptive  optics,  that  of  jitter  and  that  of  partial  coherence  of  the 
source.  We  will  now  assemble  these  results  and  propose  a  method  of 
solution. 


In  the  treatment  of  complex  problems  such  as  the  one  considered 
in  this  report,  it  is  important  to  nondimensionalize  the  governing 
equations.  This  generally  permits  a  substantial  reduction  of  the  number 
of  parameters  on  which  the  solutions  depend;  actually,  the  number  of 
parameters  is  thus  reduced  to  its  essential  minimum.  There  results  a 
saving  in  computational  effort:  the  same  solution  being  applicable  to 
different  physical  situations.  The  solutions  of  the  nondimensionalized 
system  are  called  similarity  solutions  and  the  nondimensional  parameters, 
the  similarity  parameters. 


The  propagation  distance  z  is  normalized  by  the  length  scale 


—  n.  7" 


zA=  n, 


l  -In- 


[72] 


the  transverse  coordinate,  by  the  radius  r  of  the  source  irradiance 
profile;  the  average  amplitude,  by  the  square  root  of  the  on-axis 
source  irradiance;  and  the  second-order  amplitude  moments,  by  the 
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on-axis  source  irradiance.  Performing  the  normalization  of  the  dependent 
and  independent  variables  in  eqs.  7  to  9,  we  obtain  the  following  set  of 
nondimensional  partial  differential  equations: 


>  -  b R(r,)Vp'  <  A  >  =  -  ;  ~  <  A 


{dr)  Tj  -  f  dp!  2  '•  (rj  -  f) 

(~k  +  w  ^ )  <  “*  >  “  (T^V)  f  ^  c  aa*  >  -£ Real  [ 1  <  aa* 


[73] 


2(^~aa,f~]>  +  Real  [  A(t))  ]  <  A  >  <  A  >*  +  2b  Real  [  R(rj)  ]  V„  <  A  >  ■  V„  <  A  >*, 

[74] 


/  3  p  d  \  jb  r-,2 

(  —  +  - )  <  aa  >  -  *-  Va 

Vd-rj  it  -  f  dp  /  4 

+  K(  T))  <A><A>  +  2b  /?(t))  Vp  <  A  >•  Vp  <  A  > 


Idtj  r/  -  f  dp  , 

2  <  aa  > 
<t?  -  0 


<  aa  >  -  —R(ri)  Vp  <  aa  >  +  j7'(T7)  <  aa  > 


The  coefficients  Tfn).  Cn)  and  £(n)  are  given  by 

f _  n/7  f  !/»?;»»<  1  >. 

7T77)  =  (0.15  -  0.01  j)~ -liiJ  l 

**  ^  n  n  1  J 


where 


»i ;  n  >  i 

R(b)  =  «c(»?)  +  Ra(rt)  +  «b(J?)  +  *,(>7)  -  R2(77)  -  *,(>/)  +  Rf(v) . 
K(V)  =  *„(»?)  +  AT,(»7)  -  *2(>?)  -  +  A-,(77) , 


*c(»j;f.P)  =  4  (1) 


Rn(rr,  f,j)  = 


2  f  —  >7 
J 


2  f-ij  ’ 


f  f  (2XZ  -  Y2)q  -  XY 
r-ij)\  fZQ(X  +  Yi)  +  Z172) 


Kb(rt  ;  f,  jj,,  B)  = 


jB 

2  (f- 
(Y?j  +  2X) 


[75] 


[76] 

[77] 

[78] 

[79] 

[80] 

[81] 


_  f—  n  t  Y  2(f  —  V)  t  u 

Q(X  +  Y?j  +  Zt)2)  f  fZQ  Qf  2fZ 

fY  (f- r;)  2X  Y  f-  »T|  VI  [82] 

LQ  f  fQ  2fZ  f  J  Wj  ' 
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with 


and 


with 


Hl  (f  ~~  ^ 

2  P  ' 

_  _3  yih_  (f-  n) 
2  f  f 


Z  =ji?,/f. 

Q  =  Y2  -  4XZ  . 

U  =  >n  [<X  +  Y r,  +  Zn2)/X]. 


K,(t);  f.  T)„ ,  r),  >  =  3.24 
R,(t).  f.  t)„,  yj |  )  =  2.43 


(f  -  x)-  G(y, ) 
f, 7/h  (f  -  tj »  (f  +  7)  -  2y)’ 


_ (f  -  x)-’  H(y, ) 

f 1,6  (f  -  tj)  (f  •+■  7)  -  2 y ) 


I,  =  T)„  +  0.199^-^'  t,,; 

I*  =  f,  =  rj„  +  0.335  <f  ~  n, ;  l,  =  i?„ 

y.  =  f;  y2  =  tj;  y.i  =  y*  =  y- 
v  =  _  :  L  II  -  y>  <f  ±  v  -  2y>  . 

Y'  1  2(7}  -  y)  (f  -  X)- 


The  nondimensional  independent  variables  are 

7)  =  z/zA ;  />  =  r/r„  . 

The  beam  is  characterized  by  the  similarity  parameters 

f  =  F/za  ;  b  =  z.A  / kr,:,  . 


[83] 

[84] 


[85a, b] 

f  86a, b] 


ii 
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the  turbulent  medium,  by  the  parameter 

Vo  =  0.123  ktf,/zA, 

the  adaptive  optics,  by  the  parameter 

V,  =  krf/zA 

the  noise  of  the  active  elements,  by  the  parameter 

B  =  kzA<rg,  [89] 

the  jitter  of  the  transmitting  optics,  by  the  parameter 

J=kzAff2,  [90] 

and  the  partial  coherence  of  the  source,  by  the  parameter 

P  =  2zA/kr?.  [91] 


[87] 

[88] 


The  differential  eqs.  73  to  75  together  with  the  constitutive 
eqs.  76  to  84  for  the  closure  coefficients  constitute  the  propagation 
model  of  this  report.  These  equations  are  sufficient  to  solve  for  the 
average  irradiance,  eq.  13,  and  the  irradiance  variance,  eq.  14,  of 
laser  beams  propagating  in  turbulence  and  corrected  by  adaptive  optics; 
they  also  include  modeling  of  the  effects  of  jitter  and  partial  coher¬ 
ence.  Solutions  can  be  computed  after  the  initial  profile  and  the 
various  similarity  parameters  defined  by  eqs.  86  to  91  have  been  specified. 

General  solutions  of  eqs.  73  to  75  can  be  worked  out  in  terms  of 
Green's  functions  but  this  leads  to  complicated  multiple  integrals 
that  are  troublesome  and  lengthy  to  evaluate  for  specific  applications. 

We  find  it  much  more  convenient  to  solve  the  finite-difference  version 
of  eqs. 73  to  75.  This  approach  is  straightforward  and  involves  no 
particular  difficulties  except  perhaps  in  the  near  vicinity  of  n  =  f. 
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For  originally  Gaussian  beams,  we  find  that  the  irradiance  pro¬ 
file  is  very  nearly  self-similar.  This  interesting  property  suggests 
an  algorithm  that  minimizes  the  computational  difficulties  related  to 
the  change  of  the  transverse  scale  with  the  propagation  distance.  Indeed, 
for  beams  focused  at  n  =  f,  the  diameter  decreases  with  n  until  it 
reaches  a  minimum  value  at  some  n  <  f.  If  correction  is  effective,  the 
minimum  diameter  can  be  very  small  relative  to  the  characteristic 
radius  rQ.  Hence,  solving  eqs.  73  to  75  on  a  homogeneous  grid  requires 
that  the  latter  be  chosen  fine  enough  to  resolve  the  sharpest  profile 
and  the  method  becomes  inefficient  for  the  greatest  part  of  the 
propagation  domain.  This  is  further  complicated  by  the  fact  that  the 
value  of  the  minimum  diameter  is  not  known  a  priori  except  for  the 
extremum  diffractional  limit  which  is  generally  too  small.  However, 
the  property  of  near  self-similarity  provides  a  convenient  way  around 
this  problem  as  outlined  by  the  following  method. 

The  dependent  variables  are  written  as  follows: 


<  A  >  =  F,(ij)  =  p M*)), 

[92] 

<aa*>  =  F2(i?)  P ,(»?.{). 

[93] 

<  aa  >  =  Fj(f?)  Pj(>7-  £)  • 

[94] 

If  the  solutions  were  truly  self-similar,  the  functions  would  be 
independent  of  n  and  depend  only  on  £.  However,  for  near  self-similar 
solutions,  the  functional  representation  given  by  eqs.  92  to  94  becomes 
advantageous  since  the  scale  variations  are  almost  completely  accounted 
for  by  the  functions  F^(n)  and  o(n)  of  one  independent  variable  only. 

It  is  important  to  note,  however,  that  the  eqs.  92  to  94  do  not  alter  or 
restrict  the  generality  of  the  problem,  they  are  simply  used  to  ease 
the  numerical  calculations  despite  the  apparently  more  complex  algebra. 
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Substituting  the  eqs.  92  to  94  for  <A>,  <aa*>  and  <aa>  in  the 
eqs.  73  to  75  and  assuming  circular  symmetry,  we  obtain  the  following 
equations  for  the  functions  P^(n,C): 


ap,  [Real  {*07) }  +  bfn/(f-  q)l  r  £  ap.  ,  2  P 

3i»  <t\t))  L 


] 


-b 


[/?(*/) +  j/2]  2p 
'7iP' 


=  0, 


[95] 
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[Real  {/?(*;)  }  +  bfr?  /(f  —  *;)]  f  ^  3P 
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■+  2P2  +  yV(P2 
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F2(*;) 


I2  .  ,h  Real  {*(»?)} 
F,(,)P,I  2b  o\t,)F2(v) 
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3P,  u  [Real  {*(* j) }  +  bf*?/(f-  *?)]  ["t  3P,  .  ,  D  "|  t  [R(v)  +  j/2]  „ 

_L_b - - -  ^__L+ 2  PjJ_b  V{Pj 

K(n)  Ffr>)  [971 

F,(*7)  Pl  +  ffH*?)  Fs(i7)  )  ' 


where 


^V)  =  ~~  +  bY  +  2b(f -  >;)2  J  Re(af‘ }  dr ,  [98] 

o 

1991 

Fj(*;)=  \/o\t)),  [100] 

¥}{v)=^i)exp{  -jj TTDdr}  ■  [101] 

o 


A  rapid  inspection  of  the  eqs.  95  to  97  for  the  functions  reveals 

that  the  terms  involving  differentiation  with  respect  to  5  are  all  multi¬ 
plied  by  the  similarity  parameter  b.  For  most  applications,  b  is  small 
and,  therefore,  the  resulting  functions  P^  experience  little  changes 
in  their  transverse  scale.  This  is  particularly  true  for  Gaussian 
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3P 

beams  which  satisfy  the  equation  (Stt-  +  2P.  +  1/2  V2  p  =  0)  at  n=0. 

dt,  1  i,  i. 

Consequently,  the  partial  differential  equations  95  to  97  can  be  effi¬ 
ciently  and  easily  solved  by  finite-difference  techniques  on  a  homogeneous 
grid.  The  most  important  scale  changes  occur  in  the  eqs.  98  to  101  for 
a2(n)  and  F^(n)  which  necessitate  only  three  straightforward  numerical 
integrations.  Finally,  since  the  functions  P^  are  very  nearly  self- 
similar,  the  beam  divergence  can  be  calculated  with  good  accuracy  from 
eq.  98  alone.  Inspection  of  the  latter  together  with  eq.  77  for  /?(n) 
shows  that  the  often  used  method  of  approximating  beam  divergence  by 
summing  the  squares  of  the  contributions  due  to  geometrical  optics, 
diffraction,  turbulence,  jitter  and  partial  coherence  is  verified  by 
our  model . 


5.0  RESULTS 

Our  propagation  model  was  tested  against  experimental  data  in 
Refs.  2  and  3.  Excellent  agreement  was  demonstrated  with  both  simulation 
and  atmospheric  measurements  on  widely  different  scales.  These  compari¬ 
sons  were  carried  out  for  propagation  in  turbulence  without  correction, 
jitter  or  conditions  of  partial  coherence.  In  this  report,  we  present 
calculation  results  that  simulate  the  effects  of  adaptive  optics  but 
no  data  are  yet  available  for  direct  verification.  Nevertheless,  the 
solutions  are  consistent  and  informative  and  they  show  interesting 
properties  that  should  promote  future  experimention . 

For  the  present  calculations,  we  consider  two  beams,  one  at 

10.6  ym  and  the  other  at  3.8  ym,  propagating  under  the  same  experimental 

conditions.  The  atmospheric  turbulence  is  assumed  homogeneous  with  a 

—6  —1/3 

very  strong  index  structure  constant  C  =  10  m"  '  and  an  inner  scale 

n 

lQ  =  0.01  m.  The  influence  of  the  inner  scale  on  the  solutions  is 
small  and  the  exact  value  of  a  is  not  important.  The  sources  are 
Gaussian  with  an  e-folding  irradiance  diameter  2rQ  =  0.5  m  and  the  out- 
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going  beams  are  focused  at  F  «  3  km.  The  adaptive  mirror  is  circular 
with  a  diameter  equal  to  1  m  and  contains  a  variable  number  of  active 
elements.  The  effective  area  of  the  active  elements  is  simply  taken  as 
the  total  area  of  the  mirror  divided  by  the  number  of  elements.  Finally, 
to  reduce  the  multiplicity  of  solutions,  partial  coherence  and  beam 
jitter  are  neglected.  This  has  little  consequence  for  the  purpose  of 
this  report  since  their  principal  effect  consists  in  lowering  the 
attainable  power-density  limit  with  trivial  changes  to  the  structure 
of  the  solutions. 


Figures  2  and  3  provide  examples  of  the  types  of  solutions 
calculated  with  the  present  model.  They  show  the  normalized  beam  radius 
defined  as  follows: 


CO 

p'<  I Oi,p)> 


I(i7,*>)>  dp 


[102] 


plotted  against  the  propagation  distance.  Two  curves  are  drawn:  one 
without  correction  and  the  other  with  correction.  Figure  2  pertains 
to  the  10.6-ym  wavelength  and  Fig.  3,  to  3.8  pm.  These  calculations 
show  that  the  focusing  power  at  F  =  3  km,  which  was  almost  completely 
lost  because  of  turbulence  spreading,  is  significantly  restored  by  a 
36-element  adaptive  mirror. 


To  collect  a  greater  number  of  calculation  results,  the 
Strehl  number  is  plotted  versus  the  propagation  distance  in  Figa  4  and  5 
for  10.6  um  and  3.8  um  respectively.  The  Strehl  number  is  defined  as 
the  ratio  of  the  average  on-axis  irradiance  calculated  in  turbulence  to 
the  diffraction-limited  irradiance.  A  value  of  unity  means  that  the 
correction  is  complete  and  that  the  diffraction  limit  is  attained. 

When  no  correction  is  applied,  Figs.  4  and  5  reveal  that  the  Strehl  number 
continuously  decreases  with  the  propagation  distance  up  to  the  focal 
plane  where  it  reaches  values  as  low  as  0.0082  and  0.00085  respectively. 
With  correction,  the  ratio  decreases  a  little  more  rapidly  at  the 


FIGURE  2  -  Computed  normalized  beam  radius  versus  propagation  distance 
for  a  10.6-um  Gaussian  beam.  The  other  parameters  are: 

C  =  10“^m'^'^;  lQ  =  10  mm;  rQ  =  0.25  m;  F  =  3  kmjdiameter 
or  the  deformable  mirror  D  =  1  m;  and  system  noise  =  0. 
Continuous  line:  no  correction;  dashed  line:  correction 
with  m  =  36  active  elements. 


FIGURE  3  -  Computed  normalized  beam  radius  versus  propagation  distance 
for  a  3.8-um  Gaussian  beam.  Continuous  line:  no  correction; 
dashed  line:  correction  with  m  =  36  active  elements.  Cn, 
lo<  ro*  F»  D  anc*  °0  as  in  Fig.  2. 
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beginning  but  much  more  slowly  as  z  =  F  is  approached.  If  the  number 
of  active  elements  is  sufficiently  large,  the  ratio  passes  through  a 
minimum  at  an  intermediate  propagation  distance  to  rise  again  toward  the 
focal  plane.  Ultimately,  the  diffraction  limit  is  attained  at  z  =  F 
and  little  improvement  at  shorter  distances  is  gained  by  increasing  the 
number  of  active  elements.  It  takes  more  elements  to  achieve  the  same 
diffraction-limit  condition  at  3.8  pm  than  at  10.6  pm. 

To  demonstrate  the  correction  efficiency  of  adaptive  systems, 

Figs.  6  and  7  show  the  power-density  gain  in  the  focal  plane  as  a  function 
of  the  number  of  active  elements.  The  gain  is  defined  as  the  ratio 
of  the  on-axis  focal  plane  average  irradiance  with  correction  to  the 
irradiance  calculated  without  correction.  The  gain  increases  mono- 
tonically  with  the  number  of  active  elements  and  becomes  quite  substantial 
with  relatively  few  elements.  It  is  interesting  to  note  that  the  gain 
below  the  diffraction  limit  is  of  the  same  order  at  both  wavelengths. 

The  diffraction  limit  of  about  20  dB  at  10.6  pm  is  achieved  with  60 
elements  whereas  200  elements  are  not  sufficient  to  reach  the  same 
condition  at  3.8  pm. 

Real  systems  suffer  from  noise  which  degrades  their  performance. 

To  illustrate  this  effect,  two  additional  gain  curves  are  plotted  in 
each  of  Figs.  6  and  7  for  noise  levels  =  10  urad  and  31.6  prad 
respectively.  The  parameter  cj^  is  the  standard  deviation  of  the  random 
error  in  the  wave  front  shape  returned  by  the  active  elements.  It  is 
shown  that  a  noise  below  10  prad  causes  little  degradation  but  that 
substantial  gain  reduction  occurs  at  =  31.6  urad.  These  results 
assume  the  particular  and  simple  noise  model  defined  by  eq.  34.  Although 
they  should  not  be  quantitatively  applied  to  real  systems  without  further 
verification,  they  are  indicative  of  the  magnitude  of  the  effect  of  the 
random  errors  of  the  active  elements. 


i 

i 
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FIGURE 


FIGURE 


6  -  Computed  gain  in  on-axis  focal -plane  power  density  versus  the 
number  of  active  elements  for  a  10.6-ym  Gaussian  beam, 
a:  ab  =  0;  0:  =  10  yrad;  and  A:  o.  =  31.6  yrad.  The 

horizontal  line  defines  the  diffractional  limit.  C  ,  £  ,  r  , 
F  and  D  as  in  Fig.  2.  n  o  o 


7  -  Computed  gain  in  on-axis  focal -plane  power  density  versus  the 
number  of  active  elements  for  a  3.8  ym  Gaussian  beam, 
a:  °b  s  0;  0:  %  =  io  yrad;  and  A:  av  =  31.6  yrad.  The 
horizontal  line  defines  the  diffractional  limit.  C  ,  £  ,  r  , 
F  and  D  as  in  Fig.  2.  * 
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The  effectiveness  of  the  adaptive  system  at  the  two  wavelengths, 

X  =  10.6  pm  and  3.8  pm,  is  compared  in  Figs.  8,  9  and  10  for  the  noise 
levels  =  0,  10  and  31.6  prad  respectively.  In  these  figures,  the 
on-axis  focal-plane  average  irradiance  normalized  by  the  on-axis  trans¬ 
mitter  irradiance  is  plotted  versus  the  number  of  active  elements. 

The  conditions  of  operation  are  the  same  at  both  wavelengths  so  that 
the  comparison  is  performed  in  terms  of  absolute  quantities.  In  general, 
the  turbulence-induced  beam  spread  is  greater  at  shorter  wavelengths. 

As  shown  in  Figs.  8-10,  this  remains  valid  in  the  presence  of  adaptive 
optics  up  to  a  certain  limit.  Indeed,  the  10. 6 -pm  beam  consistently 
gives  the  highest  focal-plane  power  density  until  the  10.6-pm  diffraction 
limit  is  attained.  Beyond  that  point,  increasing  the  number  of  active 
elements  no  longer  affects  the  10.6-pm  beam  but  continues  to  improve 
the  performance  of  the  shorter  3. 8- pm  wavelength  beam.  The  cross-over 
point  where  3.8  pm  becomes  more  advantageous  than  10.6  pm  is  approxi¬ 
mately  60  elements  at  =  0,  90  elements  at  =  10  prad,  and  greater 
than  200  elements  at  =  31.6  prad.  Of  course,  these  numbers  apply 
to  the  special  case  studied  in  this  report.  For  weaker  turbulence,  the 
model  predicts  that  the  cross-over  point  occurs  at  a  smaller  number  of 
active  elements. 

The  on-axis  focal  plane  scintillation  strength  is  plotted  in 
Figs.  11-13  versus  the  number  of  active  elements  for  =  0,  10  and 
31.6  prad  respectively.  The  scintillation  strength  is  defined  as  the 
ratio  aj/<I>  where  is  the  irradiance  standard  deviation  at  one  point 
and  <I>,  the  average  irradiance  at  the  same  point.  Two  curves  are 
drawn  in  each  figure  corresponding  to  the  two  wavelengths  X  =  10.6  and 
3.8  pm.  These  results  show  that  the  scintillation  remains  strong  even 
when  the  average- irradiance  beam  profile  is  near  diffraction  limited; 
a  reduction  relative  to  the  no-correction  level  of  at  best  25%  is 
observed  for  the  case  of  Figs.  11-13.  Such  residual  irradiance  fluctuations 
are  experimentally  observed,  for  example,  in  Fig.  6  of  Ref.  11. 
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FIGURE  8  -  Comparison  of  the  computed  on-axis  focal -plane  average  irra- 
diance  versus  the  number  of  active  elements  for  the  two  wave¬ 
lengths  X  =  10.6  pm  (□)  and  X  =  3.8  um  (0)  without  system 
noise,  i.e.  ab  =  0.  The  horizontal  lines  define  the  respective 

diffractional  limits.  C  ,  £  ,r  ,  F  and  D  as  in  Fig.  2. 

n  o  o  & 


FIGURE  9  -  Comparison  of  the  computed  on-axis  focal -plane  average 

irradiance  versus  the  number  of  active  elements  for  the  two 
wavelengths  X  =  10.6  pm  (□)  and  X  =  3.8  (0)  with  system 
noise  a.  =  10  prad.  The  horizontal  lines  define  the  respec¬ 
tive  diffractional  limits.  C  ,  2  ,  r  ,  F  and  D  as  in  Fig.  2. 
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FIGURE  10  -  Comparison  of  the  computed  on-axis  focal -plane  average 

irradiance  versus  the  number  of  active  elements  for  the  two 
wavelengths  X  =  10.6  pm  (□)  and  X  =  3.8  pm  (0)  with  system 
noise  =  31.6  prad.  The  horizontal  lines  defines  the 
respective  diffractional  limits.  C  ,  £  ,  r  ,  F  and  D  as  in 
Fig.  2.  °  o 

Theoretically,  a  nonzero  scintillation  level  is  not  incompatible  with  a 
diffraction- limited  beam  profile;  irradiance  fluctuations  can  indeed 
propagate  on  a  smooth  phase  front.  Moreover,  it  can  be  shown  from  the 
model  equations  that  correction  for  both  beam  spreading  and  scintillation 
would  require  the  adaptive  system  not  only  to  reproduce  the  conjugate 
of  the  phase  of  the  spherical  wave  retro-ref lected  by  the  target  glint, 
but  to  return  its  full  ’spatial'  complex  conjugate  including  the 
amplitude  modulations.  This  is  possible  with  the  nonlinear  optical 
technique  of  phase  conjugation  (e.g.  four-wave  mixing)  which  has  been 
demonstrated  in  recent  years  (Ref.  12).  However,  many  developments  are 
still  required  for  the  method  to  be  technically  applicable  to  high-power 
laser  beams . 
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Finally,  the  10.6-pm  curves  of  Figs.  8-9  and  11-12  show  that  the 
scintillation  is  minimum  at  the  number  of  elements  just  needed  to 
achieve  diffraction-limited  beam  spread.  Therefore,  the  latter  condition 
is  the  optimum  operating  condition  for  the  type  of  adaptive  systems 
studied  here.  Increasing  the  number  of  active  elements  not  only  produces 
no  further  gain  in  power  density,  but  it  worsens  the  level  of  residual 
scintillation. 


FIGURE  11  -  Comparison  of  the  computed  on-axis  focal -plane  scintillation 
strength,  Oj/<I>,  versus  the  number  of  active  elements  for 
the  two  wavelengths  X  =  10.6  pm  (□)  and  X  =  3.8  pm  (0) 
without  system  noise,  i.e.  oK  =  0.  C  ,  i  ,  r  ,  F  and  D 
as  m  Fig.  2. 
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1  10  100  200 


Number  of  Active  Elements 

FIGURE  12  -  Comparison  of  the  computed  on-axis  focal-plane  scintillation 

strength,  Oj/<I>,  versus  the  number  of  active  elements  for 

the  two  wavelengths  A  =  10.6  ym  (□)  and  A  =  3.8  um  (0)  with 

system  noise  a,  =  10  yrad.  C  ,  &  ,  r  ,  F  and  D  as  in  Fig.  2. 
d  n  o’  o’ 


1  10  100  20C 


Number  of  Active  Elements 

FIGURE  13  -  Comparison  of  the  computed  on-axis  focal-plane  scintillation 
strength,  Oj/<I>,  versus  the  number  of  active  elements  for 
the  two  wavelengths  A  =  10.6  um  (□)  and  A  =  3.8  um  (0)  with 
system  noise  =  31,6  yrad.  C  ,  l  ,  r  ,  F  and  D  as  in 
Fig.  2. 
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6.0  CONCLUSION 

The  calculation  results  of  the  preceding  section  show  that  our 
turbulence  propagation  model  can  efficiently  predict  the  average  irra- 
diance  and  the  irradiance  variance  profiles  of  adaptively  corrected 
laser  beams  traveling  in  atmospheric  turbulence.  The  solutions  are 
computed  by  a  straightforward  finite-difference  technique  and  are 
uniformly  valid  at  arbitrary  scintillation  levels. 

The  cases  studied  demonstrate  the  potential  of  adaptive  optics 
to  correct  for  turbulence  beam  spreading.  Gain  in  focal -plane  power 
density  of  the  order  of  20  dB  relative  to  uncorrected  propagation  appears 
possible  with  a  practical  number  of  active  elements  even  in  strong 
turbulence  and  in  the  presence  of  system  noise.  Gains  of  this  magnitude 
are  certainly  advantageous  whatever  the  application.  The  model  should 
thus  be  very  useful  for  the  design  of  practical  systems  and  for  scaling 
test  results  to  expected  operational  conditions. 

One  aspect  of  adaptive  optics  not  fully  accounted  for  in  the 
proposed  model  is  phase  retrieval.  Here,  we  have  assumed,  as  shown  by 
eq.  18,  that  the  phase  front  angle  returned  by  the  adaptive  system  is 
the  spatial  average,  over  the  active  element  surface  area,  of  a  quantity 
which  is  linearly  related  to  the  true  angle  of  the  conjugate  phase  front 
of  the  spherical  wave  reflected  by  the  target  glint.  In  real  applica¬ 
tions,  the  relation  can  be  more  complex  than  eq.  18  and  can  also  be 
affected  by  scintillation  and  target  speckles.  Thus,  eq.  18  will  likely 
need  more  investigation  if  it  is  applied  to  specific  adaptive  systems. 
However,  for  the  time  being,  it  remains  a  satisfactory  approximation  to 
work  out  estimates  of  the  corrective  effects  of  adaptive  optics. 
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